Period Doublings in Coupled Parametrically Forced Damped 

Pendulums 



On 
OS 



6 

o 



Sang-Yoon Kim * 

School of Physics, Georgia Institute of Technology, Atlanta, Georgia 30332-0430 

Kijin Lee 

Department of Physics, Kangwon National University, Chunchon, Kangwon-Do 200-701, Korea 



Abstract 

> 

q ' We study period doublings in N (N = 2,3,4,...) coupled parametrically 

O" 

forced damped pendulums by varying A (the amplitude of the external driv- 

O ' 
ON 



ing force) and c (the strength of coupling). With increasing A, the stationary 



point undergoes multiple period-doubling transitions to chaos. We first inves- 



tigate the two-coupled case with N = 2. For each period-doubling transition 
to chaos, the critical set consists of an infinity of critical line segments and the 
zero-coupling critical point lying on the line A = A* in the A — c plane, where 



^ ' A* is the ith transition point for the uncoupled case. We find three kinds of 



critical behaviors, depending on the position on the critical set. They are the 

same as those for the coupled one-dimensional maps. Finally, the results of 

the N = 2 case are extended to many-coupled cases with N > 3, in which the 

critical behaviors depend on the range of coupling. 
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I. INTRODUCTION 



In recent years, much attention has been paid to coupled nonlinear oscillators. Such 
coupled oscillators are used to model many physical, chemical, and biological systems such 
as coupled p-n junctions Josephson-j unction arrays 0, the charge-density waves 0, 
chemical-reaction systems 0J], and biological-oscillation systems ||. They are known to 
exhibit period-doubling bifurcations (PDB's), saddle-node bifurcations, Hopf bifurcations, 
chaos as well as pattern formation. 

The coupled oscillators investigated in this paper are coupled parametrically forced 
damped pendulums (PFDP's). For a single damped pendulum, vertical oscillation of its 
support leads to a time-periodic variation of its natural frequency, and hence it is called a 
PFDP f||7j]. This simple PFDP shows richness in its dynamical behavior [§-10]. One of the 
interesting behaviors is that with increasing the amplitude A of the vertical oscillation, the 
stationary point undergoes multiple period-doubling transitions to chaos, which have been 
found in our recent work ||11|| . 

Here we study the critical behaviors of PDB's in N (N = 2,3,4, . . .) coupled PFDP's. 
The "coupling effect" of the nature, strength, and range of coupling on the critical behaviors 
are particularly investigated. The coupled PFDP's exhibit multiple period-doubling transi- 
tions to chaos. This is in contrast to the case of the coupled one-dimensional (ID) maps, 



in which only single period-doubling transition to chaos occurs Hl2| , pT3|1 . We first consider 
the simplest coupled case with N = 2. For each period-doubling transition to chaos, the 
critical set (set of the critical points) is composed of the zero-coupling critical point and an 
infinity of critical line segments lying on the line A = A* in the A — c plane, where A is 
the amplitude of the external driving force, c a coupling parameter, and A* the ith period- 
doubling transition point for the uncoupled case. It is found that there exist three kinds of 
critical behaviors, depending on the position on the critical set. These critical behaviors are 
the same as those of the coupled one-dimensional (ID) maps found by one of us and Kook 
fT3|j . The results for the two-coupled case are also extended to many-coupled cases with 
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N > 3. It is found that the critical behaviors for the many-coupled cases vary depending on 
whether or not the coupling is global. In the extreme long-range case of global coupling, in 
which each PFDP is coupled to all the other ones with equal coupling strength, the critical 
behaviors are the same as those for the two-coupled case, irrespectively of N. However, for 
any other nonglobal-coupling cases of shorter-range couplings, an important change occurs 
in the stability diagram of 2 n -periodic (n = 0,1,2,.. .) orbits in the A — c plane, and conse- 
quently the structure of the critical set becomes different from that for the global-coupling 
case. To the best of our knowledge, so far only the authors of the paper attempted 
to study the critical behaviors of PDB's in coupled oscillators. However, only the critical 
behaviors near the zero-coupling critical point were considered, because the existence of an 
infinity of additional critical line segments in the coupled ID maps was not known at that 
time. Moreover, any explicit numerical values of the scaling factors for the zero-coupling 
case were not obtained ||15| . 



This paper is organized as follows. We first introduce two coupled PFDP's and discuss 
their general properties in Sec. [111 Stability of periodic orbits, bifurcations, and Lyapunov 
exponents are also discussed there. In Sec. JTJ, the critical behaviors of PDB's for the two- 
coupled case are studied by varying two parameters A and c. We also extend the results of 
the two-coupled case to many coupled PFDP's in Sec. |V[ Finally, a summary is given in 
Sec. 0. 



II. STABILITY OF PERIODIC ORBITS, BIFURCATIONS, AND LYAPUNOV 
EXPONENTS IN TWO COUPLED PFDP'S 

In this section, we first discuss stability of period orbits in the Poincare map of the two 
coupled PFDP's, using the Floquet theory. Bifurcations associated with the stability and 
Lyapunov exponents are then discussed. 

Consider a system consisting of two identical PFDP's coupled symmetrically: 

x x = f{x u x u t) +g(x 1 ,x 2 ), (la) 
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x 2 = f(x 2 ,x 2 ,t) +g{x 2 ,x 1 ). (lb) 

where 

f(x, x, t) = —2-n^x — 27t(cg>q — A cos 2nt) sin 2itx (2) 

and g(x\,x 2 ) is a coupling function, obeying the condition g(x,x) = for all x. Here x is 
the angular position, 7 the damping coefficient, ujq the natural frequency of the pendulum, 
A the amplitude of the external driving force of period one, and we consider the coupling 
function g(xi,x 2 ) of the form, 

(• 

g(xi, x 2 ) = - [u(x 2 ) - u(xi)\ (3) 

where u(x) is a function of one variable, and c a coupling parameter. 

The two second-order ordinary differential equations (|l|) are reduced to four first-order 
ordinary differential equations: 

±1 = 2/1, (4a) 

yi = f(x 1 ,y 1 ,t) + g(x u x 2 ), (4b) 

x 2 = y 2} (4c) 

V2 = f(x 2 ,y 2 ,t) +g(x 2 ,x 1 ). (Ad) 

Consider an initial point z(0)[= (zi(0), z 2 (0))], where z\ = (x^yi) (i = 1,2). Then, its 
Poincare maps can be computed by sampling the points z(m) at the discrete time m, where 
m — 1,2,3,... . We call the transformation z(m) — > z(m + 1) the Poincare (time-1) map, 
and write z(m + 1) = P(z(m)). 

The four-dimensional (4D) Poincare map P has an exchange symmetry such that 

SiPSi(z) = P(z) for all z, (5) 

where S\(zx, z 2 ) = (z 2 , z\). The set of all points, which are invariant under the exchange of 
coordinates Si, forms a synchronous plane on which x\ = x 2 and y\ = y 2 . An orbit is called 
a(n) (in-phase) synchronous orbit if it lies on the synchronous plane, i.e., it satisfies 



Xi(m) = £2(771) = x*(m), yi{m) = y 2 {m) = y*(m) for all m. (6) 

Otherwise, it is called an (out-of-phase) asynchronous orbit. Here we study only the syn- 
chronous orbits. They can be easily found from the uncoupled PFDP, because the coupling 
function satisfies g(x*, x*) = 0. Note also that for the cases of these synchronous orbits, the 
4D Poincare map P also has the inversion symmetry such that 

S 2 PS 2 {z) = P(z) for all z, (7) 

where ^(z) = — z. If a synchronous orbit {z(m)} of P is invariant under S 2 , it is called a 
symmetric orbit. Otherwise, it is called an asymmetric orbit and has its "conjugate" orbits 
S 2 {z(m)}. 

We now study the stability of a synchronous periodic orbit with period q such that 
P q (z(0)) = z(0) but P j (z(0)) ^ z(0) for 1 < j < q - 1. Here P k means the fc-times 
iterated map. The linear stability of the g-periodic orbit is determined from the linearized- 
map matrix DP q (z(0)) of P q at an orbit point z(0). Using the Floquet theory ||16|| , the 



matrix DP q can be obtained by integrating the linearized differential equations for small 
perturbations as follows. 

Stability analysis of an orbit can be conveniently carried out in a set of new coordinates 
{X l ,Y 1 ,X 2 ,Y 2 ) defined by 

X, = , Yl = ^i±M, (8a) 

x 2 = feizf*), Y2 = (n_yA. (8b) 

Here the first and second pairs of coordinates Z\ and Z 2l defined by Zi = (Xj, Yj) (i = 1,2), 
correspond to the synchronous and asynchronous modes of the orbit, respectively. For 
example, for a synchronous orbit Z\ = (x*,y*) and Z 2 = (0,0), while for an asynchronous 
orbit, Z 2 7^ (0,0). Hereafter, we will call Z\ and Z 2 the synchronous and asynchronous 
modes of the orbit, respectively. 

Let Z(t) [= (Zi,Z 2 )] be a solution lying on the closed orbit corresponding to a syn- 
chronous q— periodic orbit with Z(t) = Z(i + q). In order to study the stability of the 



synchronous closed orbit, we consider an infinitesimal perturbation 5Z [= (SZi, 5Z 2 )] to the 
orbit. Note that bZ\ and 5Z 2 are the synchronous and asynchronous modes of the perturba- 
tion to the synchronous orbit, respectively. Linearizing the ordinary differential equations 
(01) (expressed in terms of the new coordinates) about the orbit, we obtain 



y 8Z 2 j 



J(t) 



\5Z 2 j 



where 



J 



Jo 



Here is the 2x2 null matrix, and 



Ji(t) 



\ 



f 

yfi(x*,x*,t) f 2 (x*,x*,t) j 



K f l (x*,x*,t)-2G(x*) f 2 (x*,x*,t) j 



where 



h{x*,x*,t) = 
f 2 (x*,x*,t) = 



df{x*,x*,t) 

dx* 
df(x*,x*,t) 

dx* 



-47r 2 (ul - A cos 2-Kt) cos lux* it) , 

-27T7, 



and 



G(x*) = 



dg(x 1 ,x 2 ) 



dx 2 



2 v ' 



(9) 



(10) 



(11) 



(12) 



(13) 
(14) 



(15) 



Xl=X2=X* 

Here the prime denotes the differentiation. 

Since the 4x4 matrix J of Eq. ( [TOD is decomposed into two 2x2 submatrices J± and 
J 2 , Eq. @ is reduced to two independent equations, 



SZi = Ji(t)5Zi for i = 1,2. 



(16) 
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That is, 5Z% (synchronous-mode perturbation) and 5Z 2 (asynchronous-mode perturbation) 
become decoupled for the case of a synchronous orbit. Note also that each J, (i = 1,2) is 
a g-periodic matrix. Let Wi(t) = (w^(t),wl (t)) be a fundamental solution matrix with 
Wi(0) = I. Here (t) and wf\t) are two independent solutions expressed in column 
vector forms, and I is the 2x2 unit matrix. Then a general solution of the g-periodic 
system has the following form 

8Zi(t) = Wi(t)6Zi(p), Wi(0) = I. (17) 

Substitution of Eq. (|17D into Eq. (|16D leads to an initial- value problem to determine Wi(t), 

Wi(t) = Ji(t)Wi{t), Wi(Q) = I. (18) 

In this system of new coordinates, the linearized-map matrix M (= DP q ) has the following 
block-diagonalized form, 



M 



I M x » 

y M 2 j 



(19) 



where each 2x2 submatrix Mj [= [i = 1,2) is calculated through integration of 

Eq. (0) over the period q. In order to determine the eigenvalues of M, it is sufficient to 
solve the eigenvalue problems for the two submatrices Mi and M2, independently. Here 
the submatrices Mi and M 2 determine the stability of the synchronous orbit against the 
synchronous-mode and asynchronous-mode perturbations, respectively. Note also that the 
first submatrix Mi is just the linearized Poincare map of the PFDP fllfl , and the coupling 
affects only the second submatrix M 2 . 

The characteristic equation of each submatrix Mj (i — 1, 2) is 

- trMj Xi + det M { = 0, (20) 

where trMj and detMj denote the trace and determinant of Mi, respectively. As shown in 



17 , detM,- is calculated from a formula 
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detMi = e fo tvJ > dt . (21) 

Substituting the trace of Jj (i.e., trJ, = — 2^7) into Eq. (pTj), we obtain 

det Mi = det M 2 = e" 2 " 79 . (22) 

Hence, both the submatrices Mi and M 2 have the same constant Jacobian determinant (less 
than unity). The eigenvalues, A^i and A ij2 , of are called the Floquet stability multi- 
pliers of the synchronous g-periodic orbit. The first (second) pair of stability multipliers 
(Ai,i, Ai j2 ) [(A 2 ,i, A 2j2 )] of Mi (M 2 ) is associated with stability against the synchronous-mode 
(asynchronous-mode) perturbation, and hence it may be called the pair of synchronous 
(asynchronous) stability multipliers. Note also that the pair of synchronous stability multi- 
pliers is just the pair of stability multipliers of the uncoupled PFDP ||11|| , and the coupling 
affects only the pair of asynchronous stability multipliers. 

Each pair of stability multipliers (A^i, Aj )2 ) (i = 1,2) lies either on the circle of radius 
e~ 7T ' yq , or on the real axis in the complex plane. The synchronous orbit is stable when it 
is stable against both the synchronous- mode and asynchronous- mode perturbations, i.e., 
the moduli of all its four stability multipliers are less than unity. We first note that all 
the stability multipliers never cross the unit circle in the complex plane, and hence Hopf 
bifurcations do not occur. Consequently, the synchronous orbit can lose its stability only 
when a multiplier decreases (increases) through —1 (1) on the real axis. 

Associate with each pair of stability multipliers (A^i, A.^) a quantity Ri, called the 
residue, 

1 + detM, - trMj 
Ki ~ 2(l + detM,) ' {Z6) 

which was introduced in |TjJ to characterize stability of periodic orbits in 2D dissipative 
maps with constant Jacobian determinants. Here the first and second residues Ri and 
R2 are associated with stability of the synchronous orbit against the synchronous-mode 
and asynchronous-mode perturbations, respectively. Hereafter, they will be called the syn- 
chronous and asynchronous residues, respectively. Note also that the synchronous residue 

8 



Ri is just the residue of the uncoupled PFDP [fy], and the coupling affects only the asyn- 
chronous residue R 2 - A synchronous periodic orbit is stable when < Ri < 1 for i = 1,2; at 
both ends of Ri = and 1, the stability multipliers A,'s are 1 and —1, respectively. When 
each residue Ri decreases through (i.e., Aj increases through 1), the periodic orbit loses 
its stability via saddle-noddle or pitchfork or transcritical bifurcation. On the other hand, 
when Ri increases through 1 (i.e., Aj decreases through —1), it becomes unstable via PDB, 
also referred to as a flip or subharmonic bifurcation. For each case of the PDB's and the 
pitchfork bifurcations (PFB's), two types of supercritical and subcritical bifurcations occur. 
(For more details on bifurcations in 2D dissipative maps, refer to Ref. ||19|| .) 

The stable region of a synchronous periodic orbit in the A — c plane is bounded by 
bifurcation lines associated with PDB's and PFB's (i.e., those curves determined by the 



equations Ri = and 1 for i = 1,2), as will be seen in Sec. |T|. When the boundary 
lines on which R\ (R2) = and 1 are crossed, the synchronous orbit loses its stability via 
synchronous (asynchronous) PFB and PDB, respectively. For each case of the synchronous 
(asynchronous) PFB and PDB, two types of supercritical and subcritical bifurcations take 
place. In the supercritical case of the synchronous (asynchronous) PFB and PDB, the 
synchronous orbit loses its stability, and gives rise to the birth of a pair of new stable 
synchronous (asynchronous) orbits with the same period and a new stable synchronous 
(asynchronous) period-doubled orbit, respectively. However, in the subcritical case of the 
synchronous (asynchronous) PFB and PDB, the synchronous orbit becomes unstable by 
absorbing a pair of unstable synchronous (asynchronous) orbits with the same period and 
an unstable synchronous (asynchronous) period-doubled orbit, respectively. 

Finally, we briefly discuss Lyapunov exponents of a synchronous orbit in the Poincare 



map P, characterizing the mean exponential rate of divergence of nearby orbits |20| . The 
synchronous and asynchronous modes of a nearby orbit are decoupled, because the linearized 
Poincare map DP at the synchronous orbit is just the block-diagonalized matirx M of 
Eq. ( |T9"D with q = 1. Therefore, the 2x2 submatrices Mi and M 2 of M determine the pairs of 
synchronous and asynchronous Lyapunov exponents (01,1, 01,2) and (02,1, 02,2), characterizing 



the average exponential rate of divergence of the synchronous and asynchronous modes of 
a nearby orbit, respectively, where a^i > 0^2 for i = 1,2. Since the two submatrices have 
the same constant Jacobian determinant of Eq. (|22|), each pair of the Lyapunov exponents 
satisfy 00 + 0^2 = — 27r7 [i = 1,2). Note also that the first pair of synchronous Lyapunov 
exponents (01,1, °i,2) is j us ^ the pair of the Lyapunov exponents of the uncoupled PFDP, and 
the coupling affects only the second pair of asynchronous Lyapunov exponents (02,1, 02,2) • 



III. CRITICAL BEHAVIORS IN TWO COUPLED PFDP'S 

In this section, by varying the two parameters A and c, we study the critical behaviors 
of the synchronous supercritical PDB's in the two coupled PFDP's ([!]) for 7 = 0.1 and 
ojq = 0.5. The two coupled PFDP's exhibit multiple period-doubling transitions to chaos, 
which is in contrast to the case of the coupled ID maps with only single period-doubling 



transition to chaos |T^|T^]. For each transition to chaos, the zero-coupling critical point and 
an infinity of critical line segments constitute the critical set in the A — c plane. There exist 
three kinds of critical behaviors, depending on the position on the critical set. These critical 
behaviors are found to be the same as those of the two coupled ID maps [|I3[|. 
We consider a linearly-coupled case in which the coupling function (|3|) is 

g(xi,x 2 ) = -(x 2 -xi). (24) 

Figure [I] shows the stability diagram of the synchronous orbits with low period q = 1,2. 
The stable region of a synchronous orbit is bounded by its PDB and PFB lines. The 
horizontal (non-horizontal) solid and short-dashed boundary lines correspond to synchronous 
(asynchronous) PDB and PFB lines, respectively. (Each bifurcation may be supercritical 
or subcritical.) Note also that the horizontal synchronous PDB or PFB lines extend to the 
(plus) infinity (c = 00). For the sake of convenience, only some parts (up to c = 7) of the 
infinitely long lines are drawn in the figure. 

We first consider the bifurcations associated with stability of the stationary point 
[(xi, yi, x 2 , 2/2) = (0,0,0,0)]. Its stable region is denoted by SP in Fig. |I[ When the non- 
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horizontal short-dashed line of the SP is crossed, the stationary point becomes unstable via 
asynchronous subcritical PFB. However, at the horizontal solid boundary line, it loses its 
stability via synchronous supercritical PDB, and a new stable synchronous orbit of period 2 
appears. The 2-periodic orbit is a symmetric orbit with respect to the inversion symmetry S2, 
and its stable region is denoted by SP2 in Fig. |T[ When the horizontal (non-horizontal) short- 
dashed boundary line of the SP2 is crossed, the symmetric 2-periodic orbit loses its stability 
via synchronous (asynchronous) supercritical PFB, and a pair of new stable synchronous 
(asynchronous) orbits with period 2 appears. Note that the new pair of synchronous or- 
bits is a conjugate pair of asymmteric orbits with respect to the ^-symmetry. Its stable 
region is denoted by ASP2 in Fig. [|. Unlike the cases of the lower-level stability regions (SP 
and SP2), the ASP2 is [/-shaped, because a parabolalike asynchronous PDB line also is a 
boundary line of the ASP2. An asynchronous supercritical PDB occurs at the parabolalike 
solid line, whereas an asynchronous subcritical PFB takes place at the non-horizontal short- 
dashed line. However, each synchronous asymmetric 2-periodic orbit becomes unstable via 
synchronous supercritical PDB when the horizontal solid line is crossed, and gives rise to 
the birth of a new synchronous asymmetric 4-periodic orbit. Here we are interested in such 
synchronous supercritical PDB's. 

Figure ^| shows the stability diagram of synchronous asymmetric orbits born by syn- 
chronous supercritical PDB's. Each synchronous asymmetric orbit of level n (period 2 n , 
n = 1,2,3,...) loses its stability at the horizontal solid line of its stable region via syn- 
chronous supercritical PDB, and gives rise to the birth of a synchronous asymmetric 
period-doubled orbit of level n + 1. Such an infinite sequence ends at a finite value of 
A* = 0.357 709 845 3, which is the first period-doubling transition point of the uncoupled 
PFDP ||11|| . Consequently, a synchronous quasiperiodic orbit, whose maximum synchronous 
Lyapunov exponent is zero (i.e., cr^i = 0), exists on the A = A\ line. 

We examine the treelike structure of the stability diagram in Fig. ^|, which consists of 
an infinite pile of [/-shape regions and rectangular-shape regions. Note that the treelike 



structure is asymptotically the same as that in the coupled ID maps JT3[. The [/-shape 
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branching is repeated at one side of each [/-shape region, including the c = line segment. 
The branching side will be referred to as the zero c side. However, the other side of each 
[/-shape region grows like a chimney without any further branchings [as an example, see the 
branch starting from the right side of the ASP2 in Fig. 0(b)]. As in the coupled ID maps 
I3| , this rule governs the asymptotic behavior of the treelike structure, even though there 



are a few exceptions for lower-level orbits. Other type of [/-shape regions without the zero 
c sides [e.g., the leftmost [/-shape region in the third-level stability region in Fig. 0(a)] may 
appear in the lower-level stability regions. However, the [/-shape branching for this kind of 
[/-shape region ends at some finite level, and then each side of the [/-shape region grows like 
a chimney without any further branchings. Consequently, an infinite number of successive 
branchings occur only for the case of the [/-shape region with the zero c side. 

A sequence of connected stability regions with increasing period is called a "period- 
doubling route" [fT5| . There are two kinds of period-doubling routes. The sequence of the 
[/-shape regions with the zero c sides converges to the zero-coupling point c = on the 
A = A\ line. It will be referred to as the U route. On the other hand, a sequence of 
rectangular regions in each chimney converges to a critical line segment on the A = A\ line. 
For examples, the rightmost one in Fig. |2](a) is the line segment joining the left end point 
q (= 0.343 687 • • •) and the right end point c r (= 0.484 777 • • •) on the A = A\ line, and the 
one in Fig. 0(b) is the infinitely-long line connecting the two end points q (= 4.407457- • ■) 
and c r (= oo) on the A = A\ line. This kind of route will be called a C route. Note that 
there are infinitely many C routes, while the U route converging to the zero-coupling critical 
point (A*,0) is unique. Hence, an infinite number of critical line segments, together with 
the zero-coupling critical point, constitute the critical set. 

We now study the critical behaviors on the critical set. First, consider the case of the U 
route ending at the zero-coupling critical point. We follow the synchronous orbits of period 
q = 2 n up to level n = 9 in the U route, and obtain a self-similar sequence of parameters 
(A„,, c n ), at which each orbit of level n has some given residues Ri and R 2 (e.g., R\ = 1 and 
i?2 = 0). Then the sequence {(A n , c n )} converges geometrically to the zero-coupling critical 
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point (A*, 0). In order to see the convergence of each of the two scalar sequences {A n } and 
{c n }, we define 

x — ^-A n _ Ac n 
AA n+1 Ac n+ i 

where AA n = A n — A n _i and Ac n = c n — c n -\- The sequences of 5 n and \i n are listed in 
Table |, and converge to their limit values, 5 (~ 4.67) and \i (~ —2.5), respectively. Hence 
the two sequences {A n } and {c n } obey one-term scaling laws asymptotically: 

AA n ~5- n , Ac n ~/i" n for large n. (26) 

As in the coupled ID maps, the value of the nonlinearity-parameter scaling factor 5 agrees 
well with the Feigenbaum constant (= 4.669- • ■) of the ID maps ||21|| . The value of the 
coupling-parameter scaling factor \x is also close to that (= —2.502 • • •) of the coupling- 
parameter scaling factor a of the coupled ID maps near the zero-coupling critical point |L3| . 



It has been also shown in |13j that the scaling factor a is just the first relevant "coupling 
eigenvalue" (CE) v\ of the zero-coupling fixed map of the renormalization transformation 
for the case of the coupled ID maps. In addition to V\ = a, the zero-coupling fixed map has 
another second relevant CE v 2 (= 2), which also affects the scaling associated with coupling 



in the coupled ID map |22 



In order to get a correction to the leading scaling (|26|) , we take into account the effect 
of the second relevant CE v 2 (— 2) on the scaling of the sequence {Ac n } and extend the 



simple one-term scaling law (|26|) to a two-term scaling law | 22| , |23 |: 

Ac n ~ Ci/ip + C 2 /i2 n for large n, (27) 
where | /xa | > and C\ and C 2 are some constants. This is a kind of multiple scaling law 



2§. Eq. (0) gives 

Ac n = siAc n+ i — S2Ac n+ 2) (28) 

where s x = /ii + /i 2 and s 2 = [i\\i>2- Then, fi\ and /i 2 are solutions of the following quadratic 
equation, 
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/i 2 - si/x + s 2 = 0. (29) 

To evaluate \i\ and /i 2 , we first obtain s\ and s 2 from Ac n 's using Eq. (^): 

s _ Ac n Ac n+1 - Ac n _iAc n+2 ^ _ Acl - Ac n+ iAc n _i ^ 
Ac^+i — Ac n Ac n+ 2 ^ c n+i ~~ Ac n Ac n+ 2 



Note that Eqs. p7|)-(|30|) hold only for large n. In fact, the values of Sj's and /Vs (i = 1, 2) 
depend on the level n. Therefore, we explicitly denote s^'s and /ij's by Sj jn 's and yUi,„'s, 
respectively. Then, each of them converges to a constant as n — > oo: 

lim Si „ = lim fe n = i = 1, 2. (31) 

Three sequences {yUi, n }, {^2,n}, and {/x 2 n / /i2,n} are shown in Table ||. The second column 
shows rapid convergence of the first scaling factor n\ %n to its limit value \i\ (~ —2.50), which 
agrees well with the first relevant CE v\ (= a). (Its convergence to a is faster than that for 
the case of the above one-term scaling law.) The second scaling factor /x 2 n also seems to 
converge slowly to its limit value fi (~ 3.1), whose accuracy is lower than that of fi\. As in 
the coupled area-preserving maps ||23|| , it seems from the third and fourth columns that the 
second scaling factor \i 2 may be expressed by a product of two relevant CE's v\ (= a) and 
V2 (= 2), 

A4 = -■ (32) 
It has been known that every scaling factor in the multiple-scaling expansion of a parameter 



is expressed by a product of the eigenvalues of a linearized renormalization operator [24]. 

We also study the coupling effect on the asynchronous residue i?2,n of the synchronous 
orbit of period 2™ near the zero-coupling critical point (A*,0). Figure ^| shows three plots 
of i£ 2 ,n(^.*j c ) versus c for n = 4,5, and 6. For c = 0, _R 2 ,n converges to a constant 
(= 1.300 59. . .), called the critical asynchronous residue, as n — > oo. However, when c is 

dR,2,n 



nonzero i?2n diverges as n — * oo, i.e., its slope S n (= 

ac 

critical point diverges as n — > oo. 



) at the zero-coupling 

(a*,0) 
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As in the coupled area-preserving maps [p3|l , the sequence {S n } also obeys a two-term 
scaling law, 

S n = D x vl + D 2 vl for large n, (33) 

where \u x \ > \u 2 \. This equation gives 

S n +2 = ri5 n+ i - r 2 S n , (34) 

where ri = v\ + z/2 and r 2 = v\i>2- As in the scaling for the coupling parameter, we first 
obtain r\ and r 2 of level n from Sn's: 

S n +lS n ~ S n + 2 S n -i _ $n+l ~~ SnS n +2 /Qp-\ 

~ "n+l'-'n-l '-'n '- 5 n+l'- 5 n.-l 

Then, the scaling factors z/i jra and ^2, ra of level n are given by the roots of the quadratic 
equation, 

V 2 n - r l,nVn + ?~2,n = °. (36) 

They are listed in Table |TTT| and converge to constants v\ (~ —2.503) and v 2 (~ 2) as n — > oo, 
whose accuracies are higher than those of the coupling-parameter scaling factors. Note that 
the values of v\ and v 2 agree well with those of the two relevant CE's v\ and v 2 . 

We next consider the cases of C routes, each of which converges to a critical line segment. 
Two kinds of additional critical behaviors are found at each critical line segment; the one 
critical behavior exists at both ends and the other critical behavior exists at interior points. 
In each C route, there are two kinds of self-similar sequences of parameters (A n , c n ), at which 
each synchronous orbit of level n has some given residues R\ and R 2 ; the one converges to the 
left end point of the critical line segment and the other converges to the right end point. As 
an example, consider the rightmost C route in Fig. 0(a), which converges to the critical line 
segment with two ends (A* X) ci) and (Al,c r ). We follow, in the rightmost C route, two self- 
similar sequences of parameters, one converging to the left end and the other converging to 
the right end. In both cases, the sequence {A n } converges geometrically to its accumulation 
value A\ with the ID scaling factor 5 (~ 4.67) like the case of the U route, 
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AA n ~5~ n for large n, (37) 

where AA n = A n — A n _ x . The sequences {c n } for both cases also obey the one-term scaling 
law, 

Ac n ~ jjT n for large n, (38) 



where Ac n = c n — c n _i. The sequence of the scaling factor fx n of level n is listed in Table [TV], 
and converges to its limit value /i (—2). Since the value of the coupling-parameter scaling 
factor fi is different from that (/i = a) for the zero-coupling case, the critical behavior at 
both ends differs from that at the zero-coupling critical point. We also note that the value 
of \i agrees well with that of the coupling-parameter scaling factor [y = 2) of the coupled 



ID maps near both ends of each critical line segment |T3|]. It has been also shown in [13| 
that the scaling factor v (= 2) is just the only relevant CE of a nonzero-coupling fixed map 
of the renormalization transformation for the case of the coupled ID maps. 

Figure |] shows the behavior of the asynchronous residue R2, n {A\ 1 c) of the synchronous 
orbit of period 2 n near the rightmost critical line segment in Fig. 0(a). For c = q and c r , 
i?2,n converges to a critical residue i?2 ( = 0) as n ^ oo, which is different from that for the 
zero-coupling case. The slopes S^s of i? 2 ,n at both ends obey well the one-term scaling law, 

S n ~ v n for large n. (39) 

The two sequences of the scaling factors v n of level n at both ends are listed in Table V, 
and converge to their limit values v ~ 2, which agrees well with that of the only CE [v = 2) 
of the nonzero-coupling fixed map governing the critical behavior at both ends for the case 
of the ID maps. However, for any fixed value of c inside the critical line segment, i?2,n 
converges to a critical residue i?2 (= 0.5) as n — > oo (see Fig. |). This superstable case of 
i?2 = 0.5 corresponds to the supercritical case of X% = (A^: the second critical stability 
multiplier) for the coupled ID maps |jTB[, because Eq. (|23|) of R for the case of 2D maps 
reduces to the equation of R — 0.5 * (1 — A) for the case of ID maps. We also note that as 
in the case of the coupled ID maps, there exists no scaling factor of the coupling parameter 
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inside the critical line segemnt, and hence the coupling parameter becomes an irrelevant 
one at interior critical points. Thus, the critical behavior inside the critical line segment 
becomes the same as that of the uncoupled PFDP (i.e., that of the ID map), which will 
be discussed in more details below. This kind of ID-like critical behavior was found to be 
governed by another nonzero-coupling fixed map with no relevant CE for the case of the 
coupled ID maps [ |13|] . 

There exists a synchronous quasiperiodic orbit on the A = A\ line. As mentioned in 
Sec. HI, its synchronous Lyapunov exponents are the same as the Lyapunov exponents of 
the uncoupled PFDP, i.e., ct^i = and = — 0.2tt. The coupling affects only the second 
pair of asynchronous Lyapunov exponents (cr 2i i, 02,2), characterizing the mean exponential 
rate of divergence of the asynchronous mode of a nearby orbit. The maximum asynchronous 
Lyapunov exponent o^i near the rightmost critical line segment in Fig. 0(a) is shown in 
Fig. |5[ Inside the critical line segment (q < c < c r ), the synchronous quasiperiodic orbit on 
the synchronous plane becomes a synchronous attractor with cr 2) i < 0. Since the dynamics 
on the synchronous attractor is the same as that of the uncoupled PFDP, the critical maps 
at interior points exhibit essentially ID-like critical behaviors, because the critical behavior 
of the uncoupled PFDP is the same as that of the ID maps However, as the cou- 

pling parameter c passes through q and c r , the maximum asynchronous Lyapunov exponent 
(72,1 of the synchronous quasiperiodic orbit increases from zero. Hence, the synchronous 
quasiperiodic orbit becomes unstable and ceases to be an attractor outside the critical line 
segment. Consequently, the system of the two coupled PFDP's is asymptotically attracted 
to another synchronous or asynchronous attractor outside the critical line. For example, the 
asymptotic state for c = 0.343 68 (< q) becomes an asynchronous attractor of period 512, 
while that for c = 0.484 79 (> c r ) becomes a synchronous rotational attractor of period 1. 

We also study the critical scaling behaviors of the maximum asynchronous Lyapunov 
exponent 02,1 near both ends of the rightmost critical line segment in Fig. 0(a). As shown in 
Fig. H 02,1 varies linearly with respect to c near both ends, i.e., <j 2 ,i ~ e, e = c — c* (c* = q 
or c,.). The critical exponent of 02,1 near both ends can be also obtained from the only CE 
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v = 2 of the nonzero-coupling fixed map governing the critical behavior near both ends for 
the case of the coupled ID maps. Consider a system with nonzero e (but with A = A^) near 
both ends. It is then transformed into a new one of the same form, but with a renormalized 
parameter e' under a renormalization transformation. Here the parameter e obeys a scaling 
law, 

e' = ve = 2 e. (40) 
Then the maximum asynchronous Lyapunov exponent 02,1 satisfies the homogeneity relation, 

a 2A (e') = 2a 2A (e). (41) 

This leads to the scaling relation, 

(T2,i ~ e\ (42) 

with critical exponent 

77 = ln2/liii/ = l. (43) 

As the nonlinearity parameter A is further increased from A = A*, the stationary point 
[xi = X2 = x* = 0, yi = U2 = y* = 0] undergoes a cascade of "resurrections", i.e., it will 
restabilize after it loses its stability, destabilize again, and so ad infinitum [fLl] . It was found 



m 



TTf that for luq = 0.5, its restabilizations occur through alternating synchronous sub- 
critical PDB's and PFB's, while the destabilizations take place via alternating synchronous 
supercritical PDB's and PFB's. Consequently, the two coupled PFDP's exhibit multiple 
period-doubling transitions to chaos. This is in contrast to the case of the coupled ID maps, 



in which only single period-doubling transition to chaos occurs [0.0 • 

As the first example of the multiple period-doubling transitions to chaos, we consider the 
first resurrection of the stationary point shown in Fig. 0(a). For A = A r (l) (= 3.150 509 ■ ■ •), 
a synchronous subcritical PDB occurs. Hence, the stationary point restabilizes with birth 
of a new unstable synchronous symmetric orbit of period 2 for A > A r (l). As A is increased 



from A = A r (l), the stationary point destabilizes at A = Ad(2) (= 3.224 230 • • •) via syn- 
chronous supercritical PFB, which results in the birth of a conjugate pair of synchronous 
asymmetric orbits with period 1. Fig. 0(b) shows the stability diagram of the stationary 
point and the synchronous asymmetric orbits of level n (period 2 n , n = 0,1, 2, 3, 4) near the 
c = line in the A — c plane [2S|. Each synchronous asymmetric orbit of level n becomes 



unstable at the horizontal solid line of its stable region via synchronous supercritical PDB, 
and gives rise to the birth of a synchronous asymmetric period-doubled orbit of level n + 1. 
Such an infinite sequence terminates at a finite value of A\ = 3.263 703 15 • • •, which is the 



second period-doubling transition point of the uncoupled PFDP ||11|| . Note that the treelike 
structure of the stability diagram in Fig. 0(b) is essentially the same as that in Fig. @(a). 
Hence, the critical set also consists of an infinite number of critical line segments and the 
zero-coupling critical point, as in the first period-doubling transition case. In order to study 
the critical behaviors on the critical set, we follow the synchronous asymmetric orbits up to 
level n = 7 in the U route and the rightmost C route. It is found that the critical behaviors 
are the same as those for the first period-doubling transition case. That is, there exist three 
kinds of critical behaviors at the zero-coupling critical point, both ends of each critical line 
segment and interior points. 

As the second example, we also consider the second resurrection of the stationary 
point shown in Fig. |8|(a). A synchronous subcritical PFB takes place at A = A r (2) 
(= 10.093 985 •• •). Consequently, the stationary point restabilizes with birth of a pair of 
new unstable orbits with period 1. As A is further increased, the stationary point destabi- 
lizes at A = A<i(3) (= 10.097 583 •• •) via synchronous supercritical PDB, which results in 
the birth of a new synchronous symmetric orbit with period 2. The subsequent bifurcation 
behaviors are the same as those for the first period-doubling transition case. That is, a third 
infinite sequence of synchronous supercritical PDB's follows and ends at a finite value A% 
(= 10.099 660 93- • •), which is the third period-doubling transition point of the uncoupled 
PFDP [fLl| . The third stability diagram of synchronous orbits near the c = line is shown 
in Fig. |8|(b) [p5[| . Note that its treelike structure is essentially the same as that in Fig. |2](a). 
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Hence, the critical set is composed of the zero-coupling critical point and an infinity of criti- 
cal line segments. Furthermore, the critical behaviors on the critical set are found to be the 
same as those for the first period-doubling transition case. 



In addition to the linear-coupling case (p3), we have also studied other nonlinear-coupling 



g{xi,x 2 ) = ~[x%-x?], n = 2,3. (44) 

For the first period-doubling transition case, the stability diagrams of synchronous orbits 
near the c = line for the cases of the quadratic and cubic couplings are shown in Fig. |](a) 
and H(b), respectively. Their treelike structures are essentially the same as that in Fig. 0(a). 
Hence, the zero-coupling critical point and an infinite number of critical line segments con- 
stitute the critical set for each nonlinear-coupling case. Moreover, the critical behaviors for 
these nonlinear-coupling cases are also found to be the same as those for the linear-coupling 
case. 



IV. EXTENSION TO MANY COUPLED PFDP'S 

In this section we study the critical behaviors of the synchronous PDB's in iV-coupled 
(N > 3) PFDP's in which the coupling extends to the Kth [1 < K < f (^) for even (odd) 
N] neighbor(s) with equal strength. It is found that the critical behaviors depend on the 
coupling range. In the global-coupling case, in which each PFDP is coupled to all the other 
ones with equal coupling strength, the structure of the critical set and the critical behaviors 
are the same as those for the two-coupled case, independently of N. However, for any other 
nonglobal-coupling cases, the structure of the critical set becomes different from that for the 
global-coupling case, because of a significant change in the stability diagram. 

Consider N symmetrically coupled PFDP's with a periodic boundary condition, 

Xm f{Xmi Xmi t) ~\~ (}{x m ^ X m _|_i, . . . , X m _i), Tfl 1,2,..., N. (45) 
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Here the periodic boundary condition imposes x m (t) = x m+ N(t) for all m, the function 
f(x, x, t) is given in Eq. (fj), and g(x 1: . . . , x^) is a coupling function, obeying the condition 



g(x, . . . , x) = for all x. 



A general form of coupling for odd iV (N > 3) is given by 



(46) 



K 



g(xi, ...,x N ) 



2K 



■ l=-K 



1 



K 



2K 
K = 1 



- ^ u ( x i+i) ~ u ( x i) 



l=-K 

N- 1 

? • ■ • J o 



(47) 



where c is a coupling parameter and u is a function of one variable. Here the coupling 
extends to the i^th neighbors with equal coupling strength, and the function g satisfies the 
condition (fK|). The extreme long-range interaction for K = is called a global coupling, 
for which the coupling function g becomes 

„ N 



g{x l: . . . , x N ) = — [u(x m ) - u(xi)} 

iV m=l 



N 



N 



J2u(x m ) - u(xi) 



m=l 



(48) 



This is a kind of mean-field coupling, in which each element is coupled to all the other 



elements with equal coupling strength. All the other couplings with K < 



N-l 



'e.g., nearest- 



neighbor coupling with K = 1) will be referred to as non-global couplings. The K = 1 case 
for N = 3 corresponds to both the global coupling and the nearest-neighbor coupling. 

We next consider the case of even iV (N > 2). The form of coupling of Eq. ( |4"TD holds 



for the cases of non-global couplings with K 



1,... 



N-2 
' 2 



(N > 4). The global coupling 



for K = y (N > 2) also has the form of Eq. (fi8"|), but it cannot have the form of Eq. 
because there exists only one farthest neighbor for K = y, unlike the case of odd N. The 
K = 1 case for N = 2 also corresponds to the nearest-neighbor coupling as well as to the 
global coupling, like the N = 3 case. 
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The stability analysis of an orbit in many coupled PFDP's is conveniently carried out 
by Fourier-transforming with respect to the discrete space {m} |[26|| . Consider an orbit 
{x m (t) ; m = 1,...,N} of the N coupled PFDP's (|45|). The discrete spatial Fourier 
transform of the orbit is: 

1 N 

iV m=l 

j = 0,l,...,N-l. (49) 

The Fourier transform £j(t) satisfies £j(t) = (,N-j(t) (* denotes complex conjugate), and the 
wavelength of a mode with index j is y for j < y and for j > ~. 

To determine the stability of a synchronous g-periodic orbit [xi(t) — ■ ■ ■ — xj^(t) = x*(t) 
for all t and x*(t) = x*(t + 5)], we consider an infinitesimal perturbation {5x m (t)} to the 
synchronous orbit, i.e., x m (t) = x*(t) + 8x m (t) for m = 1, . . . , N. Linearizing the iV-coupled 
PFDP's ( f45|) at the synchronous orbit, we obtain: 



df{x*,x*,t) s df(x*,x*,t) x . 

0X m + CX m 

ax ox 

N 

+Y, G i( x *) Sxi+m-u (50) 



1=1 

where 



9flf(xi, ...,a;jv) 
Gjfa;) = - 



(51) 



Hereafter the functions G/'s will be called "reduced" coupling functions of g(xi, . . . ,xn). 
Let be the Fourier transform of &E m (t), i.e., 

1 w 

<^ = ^[^(f)] = -Y,e- 2mm ^ N 6x m , 

^ m=l 

j=0,l,...,iV-l. (52) 

Here <5£o is the synchronous- mode perturbation, and all the other 5£/s with nonzero indices j 

are the asynchronous-mode perturbations. Then the Fourier transform of Eq. ( |5TD becomes: 

= df(x*,x*,t) . + g/(x%xV) 
J Si;* 5 dx* 

N 

+ £ CC^Je^P- 1 ^]^, j = 0, 1, . . . , N - 1. (53) 
z=i 
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Note that all the modes 5£/s become decoupled for the synchronous orbit. 
The equation ([53]) can also be put into the following form: 



where 



V 

I 



hit) 



, j = 0,l,...,N-l, 



(54) 



df(x* ,x*,t) 
dx 



N 



(55) 



i=i 



Note that each Lj is a g-periodic matrix, i.e., Lj(t) = Lj(t + q). Let $j(t) = (<pf\t),(j)y'(t)) 



be a fundamental solution matrix with $y(0) = /. Here <pf\t) and ^'{t) are two indepen- 
dent solutions expressed in column vector forms, and / is the 2x2 unit matrix. Then a 
general solution of the g-periodic system has the following form 



«.( 2 )/ 



$,(t) 



<%(0) 
j = 0,l,...,iV - I. 



(56) 



Substitution of Eq. flSB]) into Eq. fl54|) leads to an initial- value problem to determine $j{t), 



(57) 



Each 2x2 matrix tyj [= $>j(q)], which is obtained through integration of Eq. ( p7| ) over the 
period q, determines the stability of the q-periodic synchronous orbit against the jth-mode 
perturbation. 

The characteristic equation of each matrix tyj (j = 0, 1, . . . , N — 1) is 



Xj-ti iffjXj + detVj = 0, 



(58) 



where tr\l/, and det\l/o denote the trace and determinant of respectively. As shown in 



17 , det tyj is given by 



det = e^ trLjdt = e" 27r7<? . 



(59) 
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Hence, all the matrices \^-'s have the same constant Jacobian determinant (less than unity). 
The eigenvalues, A 3 - ( i and A 5) 2, of are called the Floquet stability multipliers, which 
are associated with the stability of the synchronous g-periodic orbit against the jth-mode 
perturbation. Since the j = case corresponds to the synchronous mode, the first pair of 
stability multipliers (A 0j i, A ,2) is called the pair of synchronous stability multipliers. On the 
other hand, all the other pairs of stability multipliers are called the pairs of asynchronous 
stability multiplies, because all the other cases of j ^ correspond to asynchronous modes. 
Like the two-coupled case [see Eq. fl23|)], we also associate with a pair of stability multipliers 
Aj ; i and Aj^ a residue Rj, 

Rj = \ T ; \ J = 0, 1, . . . , N - 1. 60 

3 2(l + dettfj) 

Here the first one Rq is associated with the stability against the synchronous-mode pertur- 
bation, and hence it may be called the synchronous residue. On the other hand, all the other 
ones Rj (j ^ 0) are called the asynchronous residues, because they are associated with the 
stability against the asynchronous-mode perturbations. 

It follows from the condition fl46|) that the reduced coupling functions satisfy 

N 

£G z (x)=0. (61) 
i=i 



Hence the matrix for j = becomes 



io(«) 



' 1 N 

df(x*,x*,t) df(x*,x*,t) 
V dx* dx* J 



(62) 



This is just the linearized Poincare map of the uncoupled PFDP |TT] . Hence the synchronous 
residue i?o becomes the same as the residue of the uncoupled PFDP, i.e., it depends only on 
the nonlinearity parameter A. While there is no coupling effect on Rq, the coupling affects 
all the other asynchronous residues Rj (j ^ 0). 

In case of the global coupling of Eq. (|48|) , the reduced coupling functions become: 

, [l-N)G(x) for/ = 1, 
Gi(x) = {' (63) 

G{x) for / ^ 1, 
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where G(x) = jfu'(x). Substituting GVs into the second term of the (2,1) entry of the 
matrix Lj(t), we have: 



N 



Y,Gi(x)e 27Ti{l ~ 1)j/N 



1=1 



for j = 0, 
— cu'(x) for j ^ 0. 

Hence all the asynchronous residues Rj (j ^ 0) become the same, i.e., R\ 



(64) 



• • — Rn-i- 

Consequently, like the two-coupled case, there exist only two independent residues Ro and 
Ri, the values of which are also independent of N. 

We next consider the non-global coupling of the form ([IT]) and define 



G(x) = 



2K + 1 



u'(x), 



(65) 



where 1 < K < 



N-2 /A^-3^ 



for even (odd) N larger than 3. Then we have 



Gi(x) = 



-2KG(x) for / = 1, 

G(x) for 2 < / < 1 + K or 

for + 1 - K < I < N, 
otherwise. 



(66) 



Substituting the reduced coupling functions into the matrix Lj(t) of Eq. (|55D , the second 
term of the (2, 1) entry of Lj(t) becomes: 

N 

E 

i=i 

where 



Y^Gi{x)e 2lTi(l ~ 1)j/N = -S N (K,j)cu'(x), 



(67) 



Sn(KJ) 



4 



2K + 



K 



sin 



njk sin(2K + 1) 



El 

N 



N 



k =i - (2K + l) S m N 



(68) 



Hence, unlike the global-coupling case, all the asynchronous residues vary depending on the 
coupling range K as well as on the mode number j. Since Sisi{K,j) = Sn(K, N — j), the 
residues satisfy 



Rj — Rn-ii j = o,i,..., at - i. 



(69) 



25 



Thus it is sufficient to consider only the case of < j < — (2^) f° r even (odd) N. 
Comparing the expression in Eq. (|67| ) with that in Eq. (|64l) for j / 0, one can easily see 
that they are the same except for the factor S^(K,j). Consequently, making a change of 
the coupling parameter c — > s n (k j) > ^ ne res idue Rj for the non-global coupling case of range 
K becomes the same as that for the global-coupling case. 

Each pair of stability multipliers (Xj i, Aj 2) {j — 0, 1, . . . , N — 1) lies either on the circle 
of radius e _7r79 , or on the real axis in the complex plane. The synchronous orbit is stable 
against the jth-mode perturbation when < Rj < 1 (i.e., the pair of stability multipliers 
(Xj,i, Xj,2) lies inside the unit circle in the complex plane). A PDB (PFB) occurs when the 
residue Rj increases (decreases) through 1 (0) [i.e., a stability multiplier decreases (increases) 
through —1 (1)]. We also note that a(n) synchronous (asynchronous) bifurcation takes place 
for j = (j 7^ 0). For more details on bifurcatios, refer to Sec. ||. 

When the synchronous residue Rq of a synchronous periodic orbit increases through 1, 
the synchronous orbit loses its stability via synchronous PDB, giving rise to the birth of a 
new synchronous period-doubled orbit. Here we are interested in such synchronous PDB's. 
Thus, for each mode with nonzero index j we consider a region in the A — c plane, in which 
the synchronous orbit is stable against the perturbations of both modes with indices and 
j. This stable region is bounded by four bifurcation curves determined by the equations 
R = 0, 1 and Rj = 0, 1, and it will be denoted by Un- 

For the case of global coupling, those stable regions coincide, irrespectively of N and j, 
because all the asynchronous residues Rj's (j 7^ 0) are the same, independently of N. The 
stable region for this global-coupling case will be denoted by Uq- Note that Uq itself is just 
the stability region of the synchronous orbit, irrespectively of N, because the synchronous 
orbit is stable against the perturbations of all synchronous and asynchronous modes in the 
region JJq. Thus the stability diagram of synchronous orbits of period 2™ (n = 1,2,3,.. .) 
in the A — c plane becomes the same as that for the two-coupled case, independently of N. 
That is, the stable regions of the synchronous orbits form a "stability tree" in the parameter 
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plane [see Figs. 2(a) and 2(b)]. Consequently, the zero-coupling critical point and an infinite 
number of critical line segments constitute the critical set. There exists one kind of critical 
behavior in the U route ending at the zero-coupling critical point, while two other kinds of 
critical behaviors exist in each C route ending at a critical line segment. The three kinds of 
critical behaviors are the same as those for the two-coupled case, independently of N. For 
more details on the critical behaviors, refer to Sec. |TJ. 

However, the stable region Un vary depending on the coupling range K and the mode 
number j for the nonglobal-coupling cases, i.e., Un = Un(K,j). To find the stability region 
of a synchronous orbit in N coupled PFDP's with a given K, one may start with the stability 
region Uq for the global-coupling case. Rescaling the coupling parameter c by a scaling factor 
s^jW]) f° r eacn nonzero 3i the stable region Uq is transformed into a stable region Un(K,j). 
Then the stability region of the synchronous orbit is given by the intersection of all such 
stable regions U^s. An important change occurs in the stability diagram of the synchronous 
orbits of period 2 n (n— 1, 2, . . .), and consequently the structure of the critical set becomes 
different from that for the global-coupling will be seen below. 

As an example, we consider the nearest-neighbor coupling case with K = 1 in four 
linearly-coupled PFDP's, in which the coupling function is given by 

g(xi,x 2 , x 3 , x 4 ) = -{x 2 + x 4 - 2xi). (70) 



Figure [T(] shows the stability regions of the synchronous 2 n -periodic (n = 1,2,3,4) orbits. 



Note that the scaling factor s ^ j. has its minimum value | at j — 2. However, for each 
synchronous orbit, 6/4(1, 2) itself cannot be the stability region, because bifurcation curves of 
different modes with nonzero indices intersect one another. We first examine the structure 
of the stability diagram in Fig. |IO|(a), starting from the left side of the stability region of 
the synchronous orbit of level 1 (n = 1). The zero c side of 6/4(1,2) including a c = line 
segment remains unchanged, whereas the other side becomes flattened by the bifurcation 
curve of the asynchronous mode with j = 1. Due to the successive flattening with increasing 
level n, a significant change in the stability diagram occurs. Of the infinite number of period- 
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doubling routes for the global-coupling case, only the U route ending at the zero-coupling 
critical point remains. Thus only the zero-coupling point is left as a critical point in the 
parameter plane. However, as shown in Fig. 0(b), the rightmost branch of the stability 
diagram, starting from the right side of the stability region of the synchronous periodic orbit 
of level 1, is the same as that for the global-coupling case except that the coupling parameter 
c is rescaled with the maximum scaling factor s ^ ^ (= 1.5) of the j = 1 mode. Hence, the 
rightmost C route ending at a critical line segment is also left. Consequently, the critical set 
for this linear-coupling case is composed of the zero-coupling critical point and one critical 
line segment. 

Consider a self-similar sequence of parameters (A n ,c n ), at which the synchronous orbits 
of period 2 n has some given residues, in the U route for the global-coupling case. Rescaling 
the coupling parameter with the minimum scaling factor 54(1,2) (= 0.75), the sequence is 
transformed into a self-similar one for the N = 4 case of nearest-neighbor coupling. Hence, 
the critical behavior near the zero-coupling critical point becomes the same as that for the 
global-coupling case. As mentioned above, the rightmost C route in Fig. |2](b) for the global- 
coupling case is also transformed into the C route in Fig. |IU|(b) for the nearest-neighbor 
coupling case by rescaling c with the maximum scaling factor S^(l, 1) (= 1.5). Hence, the 
critical behaviors at both ends and interior points of the critical line segment are the same 
as those for the global-coupling case. 

The results for the nearest-neighbor coupling case with K = 1 extends to all the other 
nonglobal-coupling cases with 1 < K < y(^y^) f° r even (odd) N. For each nonglobal- 
coupling case with K > 1, we first consider a mode with index j min for which the scaling 
factor g^Tft js becomes the smallest one and the stability region Un(K, j m m) including a c = 
line segment. Here the value of j m { n varies depending on the range K. Like the K = 1 case, 
the zero c side of Un(K, j min ) including the c = line segemnt remains unchanged, whereas 
the other side becomes flattened by the bifurcation curves of the other modes with nonzero 
indices. Thus the overall shape of the stability diagram, starting from the left zero c side 
of the stability region of the synchronous 2-periodic orbit, becomes essentially the same as 
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that for the nearest-neighbor coupling case. Consequently, only the U route ending at the 
zero-coupling critical point is left as a period-doubling route, and the critical behavior near 
the zero-coupling critical point is also the same as that for the global-coupling case. We next 
consider a mode with index j max for which the scaling factor Sn ( K ^ becomes the largest one. 
Rescaling c with the maximum scaling factor Sn ^j — y, the rightmost C route in Fig. ^|(b) 
for the global-coupling case is transformed into the C route for the nonglobal-coupling case, 
and the critical behaviors at the critical line segment are also the same as those for the 
global-coupling case. 

V. SUMMARY 

The critical behaviors of PDB's in N coupled PFDP's are investigated by varying two 
parameters A and c. As A is increased, the stationary point of the coupled PFDP's undergoes 
an infinite series of period-doubling transitions to chaos. This is in contrast to the case of 



the coupled ID maps with only single period-doubling transition to chaos |r2| , ^3| . The two- 
coupled case with N = 2 has been first studied. For each period-doubling transition to 
chaos, the zero-coupling critical point and an infinity of critical line segments constitute the 
critical set in the parameter plane. There are three kinds of critical behaviors, depending on 
the position of the critical set. They are found to be the same as those for the coupled ID 



maps []T3|| . We also extend the results of the two-coupled case to many coupled PFDP's, in 
which the critical behaviors vary depending on whether or not the coupling is global. In the 
global-coupling case, the critical behaviors are the same as those for the two-coupled case, 
independently of N. However, for any other nonglobal-coupling cases, the structure of the 
critical set becomes different from that for the global-coupling case, because of an important 
change in the stability diagram of 2 n -periodic orbits (n = 0, 1, 2, . . .). 
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TABLES 

TABLE I. In the U route, we followed a sequence of parameters (A n ,c n ) at which the pair of 
residues (i?i, n , R2,n) of the synchronous orbit of period 2 n is (1, 0). This sequence converges to the 
zero-coupling critical point (A*, 0). The scaling factors of the nonlinearity and coupling parameters 
A and c are shown in the second and third columns, respectively. 



n 5 n //„ 

2 5.286 -2.96 

3 4.692 -2.91 

4 4.665 -2.41 

5 4.666 -2.59 

6 4.667 -2.43 

7 4.670 -2.57 

8 4.665 -2.45 



TABLE II. For the case of the U route, the scaling factors /ii ift and fi2,n in the two-term scaling 
for the coupling parameter are shown in the second and third columns, respectively. A product of 
them, i s shown in the fourth column. 



n 




^2,n 




4 


-2.536 


6.87 


0.94 


5 


-2.500 


2.84 


2.20 


6 


-2.500 


2.81 


2.22 


7 


-2.504 


3.09 


2.03 
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TABLE III. The scaling factors vi >n and V2, n m the two-term scaling for the slope S n of the 
asynchronous residue i?2,n at the zero-coupling critical point are shown in the second and third 
columns, respectively. 



4 -2.599 2.783 

5 -2.511 1.923 

6 -2.503 2.004 

7 -2.503 1.998 

8 -2.503 1.999 



TABLE IV. We followed, in the rightmost C route in Fig. 2(a), two self-similar sequences 
of parameters (A n ,c n ), at which the pair of residues (Ri <n , R2,n) of the synchronous orbit with 
period 2™ is (1,0.1). They converge to both ends of the critical line segment. The scaling factors 
of the coupling paramter at the left and right ends are shown in the second and third columns, 
respectively. In both cases the scaling factors seem to converge to the same limit value fx ~ 2. 



n fl n H n 

5 1.05 3.12 

6 1.76 2.55 

7 1.85 2.26 

8 1.94 2.12 
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TABLE V. The scaling factors v n : s in the one-term scaling for the slopes S^s of the asyn- 
chronous residue i?2,n at the left and right ends of the rightmost critical line segment in Fig. 2(a) 
are shown in the second and third columns, respectively. 



n 


v n 


v n 


4 


2.156 


1.991 


5 


1.971 


2.003 


6 


2.006 


1.999 


7 


1.999 


2.000 


8 


2.000 


2.000 
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FIGURES 

FIG. 1. Stability diagram of the synchronous orbits of low period q = 1, 2 in two linearly 
coupled PFDP's. The stable regions of the stationary point, a symmetric 2-periodic orbit, and 
an asymmetric 2-periodic orbit are denoted by SP, SP2, and ASP2, respectively. The horizontal 
(non-horizontal) solid and short-dashed boundary lines correspond to synchronous (asynchronous) 
PDB and PFB lines, respectively. For other details see the text. 

FIG. 2. Stability diagram of synchronous asymmetric 2 n -periodic (n = 1, 2, 3, 4, 5) orbits of level 
n born via synchronous supercritical PDB's. ASP2 denotes the stable region of an asymmetric 
orbit of level 1, and PN also designates the stable region of an asymmetric orbit of period N 
(N= 4, 8, 16, 32). The solid and short-dashed boundary lines represent the same as those in Fig. 1. 
The stability diagram starting from the left (right) side of the ASP2 is shown in (a) [(b)]. Note its 
treelike structure. See the text for other details. 

FIG. 3. Plots of the asynchronous residue R2 jTl (A*, c) versus c near the zero-coupling critical 
point for n = 4,5, 6. 

FIG. 4. Plots of the asynchronous residue R2 jn (A*,c) versus c near the rightmost critical line 
in Fig. 2(a) for n = 5,6, 7. 

FIG. 5. Maximum asynchronous Lyapunov exponent 172,1 of the synchronous quasiperiodic 
orbit near the rightmost critical line in Fig. 2(a). The values of 0-2,1 at both ends of the rightmost 
critical line are zero, which are denoted by solid circles. 

FIG. 6. Maximum asynchronous Lyapunov exponents 02,1 of the synchronous quasiperiodic 
orbit near both (a) the left end and (b) the right end of the rightmost critical line in Fig. 2(a). 
Here e = c — c* (c* = c\ or c r ). Note that 02,1 varies linearly with respect to c near both ends. 
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FIG. 7. (a) Bifurcation diagram (plot of x* versus A) in the vicinity of the first resurrection 
of the stationary point with x* = 0; x\ = X2 = x* for a synchronous orbit. Here q = 1(2) denotes 
the period of a synchronous orbit, born via supercritical PFB (subcritical PDB). The solid and 
short-dashed lines also designate stable and unstable orbits, respectively, (b) Second stability 
diagram of synchronous orbits near the c = line. Here SP, ASP1, and PN denote the stable 
regions of the stationary point, an asymmetric orbit of period 1, and an asymmetric N-periodic 
(N= 2, 4, 8, 16) orbit, respectively. The solid and short-dashed boundary lines also represent the 
same as those in Fig. 1. For other details see the text. 

FIG. 8. (a) Bifurcation diagram (plot of x* versus A) in the vicinity of the second resurrection 
of the stationary point with x* = 0; x\ = X2 = x* for a synchronous orbit. Here q = 1(2) 
denotes the period of a synchronous orbit, born via subcritical PFB (supercritical PDB). As in 
Fig. 7(a), the solid and short-dashed lines also designate stable and unstable orbits, respectively, 
(b) Third stability diagram of synchronous orbits near the c = line. Here SP, SP2, ASP2, and 
PN (N= 4, 8, 16) denote the stable regions of the stationary point, a symmetric orbit of period 
2, an asymmetric 2-periodic orbit, and an asymmetric orbit with period N (N= 4, 8, 16) orbit, 
respectively. The solid and short-dashed boundary lines also represent the same as those in Fig. 1. 
For other details see the text. 

FIG. 9. Stability diagrams of synchronous orbits near the c = line for the cases of (a) the 
quadratic and (b) cubic couplings. Here SP2, ASP2, and PN (N= 4, 8) denote the stable regions 
of a symmetric orbit of period 2, an asymmetric 2-periodic orbit, and an asymmetric orbit with 
period N, respectively. 

FIG. 10. Stability diagram of synchronous orbits in four linearly-coupled PFDP's. Each 
stable region is bounded by its solid boundary curves. For a synchronous orbit of period q, the 
PDB (PFB) curve of the mode with index j is denoted by a symbol q^ D<yPF \ The stability diagram 
starting from the left (right) side of a 2-periodic orbit is shown in (a) [(b)]. For other details see 
the text. 
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